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Abstract 

Canine coronavirus (CCoV) exists in types I and II and infects dogs leading mainly to enteritis, though type II has 
already been associated with generalized and highly lethal infection. A CCoV-type II inactivated vaccine produced in 
A72 canine cells is available worldwide and largely used, though the molecular stability after serial passages of vaccine 
seeds is unknown. This article reports the evolution of the CCoV-II vaccine strain 1-71 in A72 cells based on partial 
S gene sequencing, showing the predominance of neutral evolution and the occurrence of four sites under purifying 
selection. Thus, cell-adapted strains of CCoV-II maybe genetically stable after serial passages in a same cell line due to 
a stable virus-host relationship. 
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Resumo 

O Coronavirus canino (CCoV) ocorre como tipos I e II e infecta caes, levando principalmente a enterite, apesar do 
tipo II ja ter sido associado a infec^ao generalizada e altamente letal. Uma vacina de CCoV-II inativada produzida 
em celulas caninas A72 e disponivel mundialmente e largamente utilizada, apesar da sua estabilidade molecular apos 
passagens seriadas de sementes vacinais ser desconhecida. Este artigo relata a evolu^ao da amostra vacinal CCoC-II 
1-71 em celulas A-72 com base em sequenciamento parcial do gene S, demonstrando predommio de evolu^ao neutra 
e a ocorrencia de quaro sitios sob sele^ao purificante. Portanto, amostras de CCoV-II adaptadas a cultivos celulares 
podem ser estaveis geneticamente apos passagens seriadas em uma mesma linhagem celular devido a existencia de uma 
rela^ao estavel virus-hospedeiro. 
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Introduction 

Canine coronavirus (CCoV), feline coronavirus 
(FCoV), and transmissible gastroenteritis virus (TGEV) are 
currently classified in the same species Alphacoronavirus-1, 


(. Nidovirales: Coronaviridae: Coronavirinae) (DE GROOT 
et al., 2013) and evolutionary studies on Alphacoronavirus-1 
have demonstrated that some geno-/serotypes of these 
viruses emerged as a result of recombination events 
between each other (PRATELLI et ah, 2004; PRATELLI, 
2006; DECARO et ah, 2009). 

CCoV presents an enveloped pleomorphic virion of 
220nm in diameter and a single-stranded positive-sense RNA 
genome of circa 29kb consisting of genes encoding the spike 
(S) envelope glycoprotein, nucleocapsid (N), membrane (M), 
and envelope (E) structural proteins and ORFs translated into 
non structural proteins, as replicase polyprotein (ORF1), 3a, 
3b, 3c, 7a, and 7b (MASTERS; PERLMAN, 2013). 
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The S protein is the main target of neutralizing antibodies 
and the most polymorphic protein in coronaviruses. It 
occurs as trimmers, with a highly variable ectodomain 
(subunit SI) and a more stable endodomain (S2 subunit) 
involved in membrane fusion (BOSCH et al., 2003; 
MASTERS, 2006). 

Estimates on the time of the most recent common 
ancestor (tMRCA) place the origin of Alphacoronaviruses 
293 million of years ago and, due to the proofreading 
activity of the NSP14 of coronaviruses, these viruses show 
a number of mutations per site per replication up to E-5, 
which is quite similar to the one found in DNA viruses 
(WERTHEIM et al.,2013). 

In spite of the large use of cell-passaged CCoV vaccine 
seeds for years on and the use of high-passage of vaccine 
strains in vaccine trials (PRATELLI, 2007), no data on 
the molecular diversity and the mechanisms of molecular 
evolution of CCoV in vitro is available. This knowledge is 
paramount for the understanding of the genetic stability 
of vaccine seeds and might also help shed some light on 
basic mechanisms of coronavirus molecular evolution. 

The aim of this investigation was to assess the clonal 
diversity of a CCoV-II vaccine strain after serial passages 
in A-72 cells based on the molecular diversity of the S2 
subunit of the spike gene and the selection regime driving 
CCoV in vitro evolution. 

Material and Methods 

CCoV strain and cell line 

CCoV 1-71 type II vaccine strain (ATCC VR 809, 
10 403 TCID 50 /mL) and canine fibroma A72 (ATCC CRL 
1542) cell line were used for this study and were kindly 
provided by Biovet Taboratories, Brazil. 

Cells were propagated in Eagle’s MEM supplemented 
with 10% fetal bovine serum (FBS) and 25mM HEPES 
(4-(2-hydroxyethyl)-l-piperazineethanesulfonic acid). 

Serial passages in A72 

For the first passage, lmL of CCV 1-71 strain diluted 
1:10 (v/v) in Eagle’s MEM was inoculated in 25 cm 2 flasks 
containing partially (-70%) confluent A72 cells and 
incubated at 37°C/ 1 hour for virus attachment. Next, the 
virus suspension was discarded and 7ml of Eagle’s MEM 
supplemented with 2% FBS was added to the flask. 

Cells were daily observed with a light microscope for 
cytopathic effect (CPE) for five days. After this observation 


period, the monolayers were frozen at -80°C regardless 
of the presence of CPE. Then, lmL was used in the next 
passage following the same protocol described above up 
to the 9 th passage 

For each passage, a 25 cm 2 with A72 cell was mock- 
infected (negative control) with Eagle’s MEM without virus 
under the same conditions described for virus inoculation. 

RT-PCR for partial S gene amplification 

For the original CCoV 1-71 and each passage (CCoV 
and negative control), total RNA was extracted with TRIzol 
Reagent " (Life Technologies) following the manufacturer’s 
instructions. cDNA synthesis was carried out with Random 
Primers'" (Life Technologies ) and M-MLV Reverse 
Transcriptase“(Life Technologies) accordingly to the 
manufacturer’s instructions. 

For partial S gene amplification, the primers 
described by Pratelli et al. (PRATELLI et al., 
2004) (S 5’TGCATTTGTGTCTCAGACTT3’/ S6 
5’CCAAGGCCATTTTACATAAG3’, nt 3492 to 4185 of 
S gene, S2 subunit in the reference strain CCoV 1-71 
GenBank JQ404409) were used with Platinum Taq 
Polymerase High Fidelity " (Life Technologies) following 
the manufacturer’s instructions. 

DNA cloning and Sequencing 

The 694 bp amplicons were purified using ExoSAP- 
IT™ (Affymetrix) and the total DNA quantified with Low 
Mass DNA Ladder " (Life Technologies) according to the 
manufacturer’s instructions. 

The purified amplicons were inserted into the plasmid 
vector pTZ57R/T'“ (InsTAclone) in an insert: vector proportion 
of 3:1 and the ligation reaction was performed at 15°C for 16 
hours, followed by the addition of 5 pL of ligation reaction to 
100 pL of competent DH5a E. coli cells grown in LB (Luria- 
Bertani) liquid medium and incubated on ice from 20 to 30 
min. After a thermal shock (42°C/2 min and ice/2 min), 900 pL 
of LB medium without antibiotics were added to the reaction 
and incubated at 37°C /150 rpm / lh30 min. 

The cells were then plated on solid LB medium 
containing 100 pg/ml ampicillin, IPTG (5-bromo-4-chloro- 
3-indolyl-(3-D-galactosidase), andX-Gal (isopropylthio-(3- 
Dgalactoside), and then incubated for 16 hours at 37°C. 

Individual white colonies, predicted to have 
inserts, were selected and subjected to a PCR reaction 
to confirm the insertion using GoTaq Green Master 
Mix'" (Promega) accordingly to the manufacturer’s 
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instructions and the plasmid-targeting primers 
M13F 5’GTTTTCCC AGT C ACG AC3 7 M13 

5' C AGGA A AC AGCTAT GAC 3’). 

A number from 9 to 13 clones for each passage and 
of the original virus confirmed as carrying the inserted 
amplicon were subjected to DNA sequencing with BigDye 
3.1“ (Applied Byosystems) and ABI-3500 Genetic Analyzer” 
(Applied Byosystems) as per manufacturers instructions, 
using primers M13F and M13R. Evaluation of the quality of 
chromatograms obtained for each DNA strand sequences was 
performed using Phred application online (http://asparagin. 
cenargen.embrapa.br/phph/) and manual checked using the 
program Finch TV (®Geospiza) to search for interpretation 
errors and discrepancies between each DNA strand 
sequenced. The final sequence of each clone was obtained 
using Cap-Contig, with PHRED score > 20, in the program 
Bioedit version 7.2.5 and submitted to BLASTn in http:// 
www.ncbi.nlm.nih.gov/BLAST for homology checking. 

Mutation analysis 

The occurrence of recombination was assessed with 
RDP, GENECONV, Bootscan, MaxiChi, Chimaera, SiScan, 
and 3Seq methods with Bonferroni correction and highest 
acceptable P-value of 0.05 using RDP 4 (336 (MARTIN et 
al., 2010). 

Sequences of clones of each passage and the original 
virus ( e.g. intra-passage analysis) were aligned using 
Clustal/W in Bioedit version 7.2.5 and compared in terms 
of nucleotide and amino acid identities. Identical sequences 
that predominated in the intra-passage comparison were 
named as Master Sequence for that passage, while those 
in lower frequencies were named Variants. 

Next, the Variants and Master Sequences for all passages 
and the original virus were aligned together including 
Alphacoronavirus homologous sequences retrieved from 
GenBank (JQ404409; JF682842; GQ477367; KC175340; 
JQ404410; EU924790; AY342160; D13096; EU924791; 
EU856361; EU856362; DQ201447; HQ462571; GQ152141; 
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NC_002306; DQ811789; DQ112226) using CLUSTAL/W 
running in Bioedit 7.2.5. 

All sequences generated in this study were deposited 
in the GenBank under the accession numbers: KP281485 
to KP281596. 

Selection regime analysis 

For each codon, estimates of the numbers of 
inferred synonymous (s) and nonsynonymous (n) 
substitutions and the numbers of sites that are estimated 
to be synonymous (S) and nonsyonymous (N) were 
produced using the joint Maximum Likelihood (ML) 
reconstructions of ancestral states under a Muse-Gaut 
model (MUSE; GAUT, 1994) of codon substitution 
and General Time Reversible model of nucleotide 
substitution. For estimating ML values, a tree topology 
was automatically computed. The test statistic dN - dS 
was used for detecting codons that have undergone 
positive selection, where dS is the number of 
synonymous substitutions per site and dN is the number 
of nonsynonymous substitutions per site. Values of p 
less than 0.05 were considered significant at a 5% level. 
These analyses were conducted using HyPhy software 
v. 2.2.1 (POND et al., 2005) implemented in MEGA6 
(TAMURA et al., 2013). 

Results 

Mutation analysis 

A total of 112 viable DNA sequences were obtained 
from passages 1 to 9 and from the original CCoV 1-71 
strain, with a number of clones per passage ranging from 
9 to 13 and only 12 nucleotide mutations found amongst 
these sequences, ranging from 0 to 3 per passage (Table 1) 
but in a passage-exclusive way. No recombination was 
found amongst these sequences. 

From these mutations, live were silent transitions 
(A3519G, T3813C, C3699T, T4155C, and A3527G); 
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three were missense transitions (A3817G, G3868A, and 
G3875A), one silent transversion (A4165C) and two 
missense transversions (G3528T andT3803G). 

Besides, a deletion (TTGATATTA/ Ile/Asp/Ile, 
nucleotides 3935 to 3943) was found in one out of 11 clones 
of passage 3. 


Selection regime analysis 

None of the 12 sites with mutations (Table 1) was 
under positive selection and only sites A3519G, T3813C, 
C3699T, and T4155C were found to be under purifying 
selection, while the remaining eight mutations were 
neutral, including the deletion. 


Table 1 - Nucleotide and putative amino acids mutations of the spike gene of CCoV stain 1-71 passaged in A72 cells. P- passage 
number (0 refers to the original virus); NC - number of clones sequenced per each passage; M/P - number of mutations 
per passage regarding the Master Sequence. Numbers in cells of each mutation represent the number of occurrences of 
that specific mutation. Positions at S gene in the reference strain 1-71 GenBank JQ404409 - Sao Paulo - 2015 


p 

NC 

M/P 

A3519G 

Arg/Arg 

A3527G 

Glu/Gly 

G3528T 

Glu/Asp 

C3699T 

Thr/Thr 

T3803G 

Phe/Cys 

T3813C 

Thr/Thr 

A3817G 

Arg/Gly 

G3868A 

Cys/Tyr 

G3875A 

Glu/Lys 

T4155C 

Asn/Asn 

A4165C 

Thr/Thr 

TTGATATTA 

3935 to 3943 
Deletion 
(lle+Asp+lle) 

0 

9 

2 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

10 

3 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

1 

0 

2 

11 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

11 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

4 

10 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

5 

13 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

6 

13 

2 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

7 

12 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

11 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

9 

12 

1 

0 

1 

0 

0 

0 
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0 

0 

0 

0 
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Figure 1 - Clones’ sequences with amino acids mutations in Alphacoronavirus-1 S2 subunit (aa 1172 to 1388 in strain 1-71 Genbank 
JQ404409) after serial passages of strain 1-71 in A72 cells 


Source: this article 


Braz. J. Vet. Res. Anim. Sci., Sao Paulo, v. 54, n. 1, p. 48-53, 2017 





















































































I 52 


Discussion 

After nine serial passages in A72 cells, the 1-71 vaccine 
strain of CCoV showed a low relative number of mutations 
when the clonal diversity of a segment of the S2 subunit 
was focused, without a passage-association pattern, which 
might indicate that this strain have reached a quite stable 
virus-cell relationship with no further mutations being 
fixed even after serial passages. 

Though the region of the S gene targeted in this 
study is responsible for membrane fusion and syncytial 
formation and has a higher mutation pressure than the 
SI region, responsible for binding to cellular receptors 
(LAI; CAVANAGH, 1997), nucleotide polymorphism may 
also be found in S2 region (GALLAGHER; BUCHMEIER, 
2001) allowing also virus entry into a variety of types of 
cells in trypsin-independent way (BELOUZARD et al., 
2009; BORUCKI et al., 2013). 

A study using a FIPV (Feline infectious peritonitis 
virus, a cat alphacoronavirus) strain passaged 100 times in 
A72 cells has shown that, from the original strain to the 8 th 
passage, amino acids substitutions occurred in SI (T599A 
and A743V), while in S2 an amino acid substitution occurred 
only at the 100 th passage (R1325E) (PHILLIPS et al., 2013). 

TGEV (Transmissible gastroenteritis virus), another 
alphacoronavirus, has been shown to have a S585A in 
SI at passage 50 in PK (Porcine Kidney) cells leading to 
attenuation of virulence (LI et al., 2010). 

As seen in Table 1, the mutations detected amongst 
the passages were not fixed, which could be due either to 
genetic drift or to a bias on the randomness of the selection 
of clones for DNA sequencing. 

Though these reports provide evidence that 
alphacoronaviruses are genetically quite stable after a high 
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